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Abstract 

We present the results of a computation of the sum of the strange and average up-down quark 
masses with overlap fermions in the quenched approximation. Since the overlap regularization 
preserves chiral symmetry at finite cutoff and volume, no additive quark mass renormalization 
is required and the results are 0{a) improved. Our simulations are performed at (3 — 6.0 
and volume ^ = 16'^ x 32, which correspond to a lattice cutoff of ~ 2 GcV and to an exten- 
sion of ~ f .4 fm. The logarithmically divergent renormalization constant has been computed 
non-perturbatively in the RI/MOM scheme. By using the if-meson mass as experimental in- 
put, we obtain (m^ + m)^^(2 GeV) = f20(7)(21) MeV, which corresponds to mf-^(2 GeV) = 
102(6)(18) MeV if continuum perturbation theory and xPT are used. By using the GMOR 

relation we also obtain {i;ij)^^{2Gey)/Nf ^ -O.Of 90(11)(33) GeV^ = -[267(5)(f 5) MeV]^. 



1 Introduction 



Quark masses are fundamental parameters of the Standard Model which cannot be measured 
directly by experiment, since quarks are confined into hadrons. If defined as effective couplings 
in the Lagrangian, their values can be determined by comparing a theoretical calculation of a 
given physical quantity (sensitive to quark masses) with the corresponding experimental value. As 
a consequence quark masses depend on the renormalization scheme and scale, as well as on the 
fundamental action. 

At present the most precise values of light-quark mass ratios (which are scheme and scale 
independent for mass independent renormalization schemes) are extracted by comparing K- and vr- 
meson mass ratios with predictions from chiral perturbation theory (xPT) A detailed analysis 
gives @ 

^ = 0.553 ±0.043, ^ = 24.4 ±1.5 (1) 

where m = (rriu ± mrf)/2. The absolute scale cannot be fixed by xPT. It can be determined 
by comparing non-perturbative lattice QCD computations or phenomenological estimates 

[0, |l^ with experimental results. 

In the last few years much effort has been devoted within the lattice community to obtain a pre- 
cise determination of the light quark masses in the quenched approximation. Major improvements 
came with the use of non-perturbative (NP) renormalization techniques for renormalizing the bare 
quark masses and with the implementation of 0{a) improved actions and operators [^, ^, 



(for a recent review see [16|). First results from unquenched simulations have also been reported 

In this paper we present the results of the first fully non-perturbative computation of (m^ ± m) 
with overlap fermions. In this fermionic regularization, flavor and chiral symmetries are preserved 
at finite lattice spacing and finite volume. As a consequence no additive quark mass renormalization 
is required and no parameters have to be fine tuned in order to compute 0{a) improved masses 
and matrix elements. To avoid uncertainties due to lattice perturbation theory, we compute the 



logarithmic divergent renormalization constant non-perturbatively in the RI/MOM scheme |17]. 
By comparing the experimental iT-meson mass with the value obtained in our simulations and 
after a careful analysis of the systematic uncertainties we find as our main result 

(m^±m)^^(2 GeV) = 120±7±21MeV (2) 

This value corresponds to 

mf^{2 GeV) = 102 ± 6 ± 18 MeV (3) 

if next-to-next-to-leading order (N^LO) continuum perturbative results and Eq. (^ are used. 

We also report results for the chiral condensate (V'V') which we compute from the Gell-Mann, 
Oakes and Renner (GMOR) relation. Our best determination is 

-^('i/iV)^H2GeV) = -0.0167 ±0.0010 ±0.0029 GeV^ (4) 

Nf 
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which corresponds to 

— (■0V')^(2GeV) = - 0.0190 ± 0.0011 ± 0.0033 GeV^ = - [267 ± 5 ± ISMeV]^ (5) 
Nf 

We have also computed the chiral condensate directly. Even if it requires a severe chiral extrapo- 
lation, the direct determination gives results consistent with Eqs. (^), (||). 

The first error in Eqs. ©-(^ represents the statistical error, obtained with the jackknife 
method. The second error represents our estimate of systematic effects. It should be noted that, 
while the systematic error due to quenching alone (defined as the quenching error in the limit 
of zero lattice spacing and infinite volume) should be the same in all fermionic regularizations, 
discretization errors do depend on the fermionic lattice action and are likely to be smaller in the 
overlap formulation, because of its good chirality properties. The very good agreement of our 
results with the current lattice world averages [16| provides further confirmation that the overlap 
formulation is a suitable regularization for large-scale phenomenological computations. 

This paper is organized as follows: in Sec. ^ we set our notation and define the renormalized 
quark masses and chiral condensate in the overlap regularization; in Sec. ^ we give details about 
the simulation and some of the results for the meson masses and matrix elements; in Sec. ^ we 
present our main results and discuss their systematic errors; Sec. ^ is devoted to our concluding 
remarks. 



2 Quark Masses and Chiral Condensate with the Overlap Action 



The QCD lattice action in the overlap regularization reads 



S 



-y 
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Up + u\, 
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—aM]D + M 
2p 



(6) 



where, in standard notation. Up \s the Wilson plaquette, go = \/6//3 is the bare coupling constant, 
ip and carry implicit color, spin and fiavor indices, and M is a diagonal matrix of bare masses 
(mi,m2, . . .) in fiavor space. D is the Neuberger-Dirac operator defined as 
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(7) 



where 

Dw = \l^.{y^. + y*^-\aVlV^ (8) 

is the Wilson-Dirac operator, < r < 1 and < p < 2r. In our calculations we used r = 1. 
and V* in Eq. (^) are the forward and backward lattice covariant derivatives, defined by 

U^{x)il){x + afi) - il){x) 



a 
1 r 



'4){x) — UUx — ajl)ilj{x — afi) 



(9) 
(10) 
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where U^{x) are the lattice gauge hnks. The fermionic operator of the overlap action satisfies the 
Ginsparg- Wilson relation |19] 

75i? + Z)75 = -D^^D (11) 
P 

which implies an exact continuous symmetry of the action in the massless limit [ 2C ] . This symmetry 
may be interpreted as a lattice form of chiral invariance at finite cutoff 

dtj) = 75^/; , 5^ = ^/i^g (12) 

where 75 is defined as 

75 = 75(1--!)) (13) 

and satisfies 

7I = 75 , 7I = 1 (14) 



The anomaly is recovered from the variation of the measure under the rotations in Eq. (12) pQ, [2l[|. 

The invariance of the action under non-singlet chiral transformations, defined including a fiavor 
group generator in Eq. (^), forbids mixing among operators of different chirality |22| and therefore: 

• No additive quark mass renormalization is required. The quark mass which enters the vector 
and axial Ward identities is the bare parameter m(a) in the action of Eq. (P). 

• Masses and matrix elements are affected only by 0(c?^ discretization errors. No fine tuned 
parameters are required to remove 0{a) effects. 

• The chiral condensate does not require subtractions of power divergent terms (in the chiral 
limit). 

The renormalized quark mass is defined as 

m(/i) = lim Zm{afJ,)m{a) (15) 

where Zm{ap) is a logarithmically divergent renormalization constant which has to be fixed, for 
a given scale p,, in a given renormalization scheme. It is worth noting that even if m(a) is the 
bare parameter which enters the fundamental action in Eq. (^), its relation to a given experimental 
result is fixed by a non-perturbative lattice QCD calculation and therefore its value is determined 
up to O(a^) terms only. 

The bare chiral condensate is defined as 

X(a) ^ hm -^(V;(0)[(1 - ^DmO)) (16) 

where m is a common mass given to the light quarks. It satisfies the integrated non-singlet chiral 
Ward identity 

^(^(0)[(1 - Y^nmO)) = mJ2{P{x)P'^m (17) 



where 



= v5i(x)75[(l--Z))V2](x) 
2p 

P\x) = i,^{x)^^[{l-^D)^l,^]{x) 

2p 



(18) 
(19) 



correspond the non-singlet pseudoscalar density with degenerate quarks (mi = m2) and its conju- 
gate. For non-zero mass the chiral condensate is still divergent and it behaves as 



^(V;(0)[(1 - Yp^)m)) = X{a) + 



(20) 



where we have taken into account that chiral symmetry forces the coefficient of the linear divergence 
to be zero. 

By writing the correlation function {P{x)P'^(0)) as a time-ordered product and by inserting a 
complete set of states in standard fashion we can also write 



Xia) 



lim 



m 



(0|P|P) 



(21) 



where Mp is the mass of the pseudoscalar state \P). If we use 

2m|(0|P|7r) I = /pM|, (22) 
where fp is the corresponding pseudoscalar decay constant, we arrive to the familiar GMOR relation 



lim 



*o 4m 



The renormalized chiral condensate is defined as 



1 



— (^V)(/^) = lini Z5(a/i)x(a) 

iV J a— >0 



(23) 



(24) 



where Zs{na) can be chosen to be the renormalization constant of the corresponding non-singlet 
scalar density which, thanks to the flavor symmetry, satisfies Zs{fJ^a) = l/Zm(/xa). In principle 
Zs{fJ.a) can be computed in perturbation theory [p3|, [2^, but uncertainties due to higher order terms 
can be avoided using non-perturbative renormalization procedures 2^, 11 1. The implementation 
of the RI/MOM technique [17| is straightforward in the overlap regularization, and it allows one to 
compute 0{a) improved renormalization constants for generic composite operators. In the following 
we will use the numerical value of Z^-^ (/xa) we have obtained in Ref . |2^] (for an alternative approach 
see Ref. pi). 
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3 Numerical Details 



We performed our simulation in quenched QCD with /3 = 6.0 and V = IG^ x 32. We used a 
sample of 54 gauge configurations, generated with the standard Wilson gluonic action of Eq. (^), 



which we retrieved from the repository at the "Gauge Connection" [28|. We computed overlap 
propagators from a local source for bare quark masses ma = 0.040,0.055,0.070,0.085,0.100 and 
p = 1.4 (and, as already mentioned, r = 1.0). For the calculation of the propagators, we used 
an optimal rational approximation to the sign of the Hermitian Wilson operator, as proposed in 
Ref. ||2^, after explicit evaluation of the contributions from the lowest eigenvectors of X^X, 
and nested multi-conjugate gradient inversions. Details of the numerical implementation will be 
presented in Ref. |26]. From the propagators, we computed in the standard manner the two-points 



correlation functionsf] 

Gss{t) = J2iS{x,t)S'iO,0)) (25) 

X 

Gpp{t) = 5](P(x,t)P^(0,0)) (26) 

X 

G^Apit) = J2i^oAo{x,t)P%0,0)) (27) 

X 

where P{x, t) and P'^{x, t) have been defined in Eqs. (|l|) and (||), 

S{x) = ^,(xm-^D)i;2]{x) (28) 
2p 

A^{x) = Vii(x)7^75[(l-^I))V'2](^) (29) 

5"^(x) is defined analogously to Eq. (|l9|) and Vq denotes the symmetric lattice derivative in the 
time direction. To improve statistics Gssit) ^-nd Gpp{t) have been symmetrized around t = T/2 
(T = Nf = 32). We estimated the errors by a jackknife procedure, blocking the data in groups 
of three configurations, and we checked that blocking in groups of different size did not produce 
relevant changes in the error estimates. In the first plot of Fig. |l| we show the ratio 

as a function of t for all simulated masses. Once p has been fitted to a constant in the time interval 
^1 — t2 = 5 — 27, a quadratic fit of the results 

(ap) = Ap + Bp{am) + Cp{am)'^ (31) 

gives^ (see the second plot in Fig. |l|) 

Ap = -0.00002(7) Bp = 1.286(3) Cp = 0.277(12) (32) 

where the quoted errors are statistical only. Note that the intercept is compatible with zero. This 

^Analogous correlation functions have been computed by the authors of Ref. A direct comparison with our 
results is not possible because we used different simulations parameters. 

^From the coefficient Bp one can derive the value of the renormalization constant Za of the "local" axial current. 
A detailed analysis will be presented in Ref. teq]. 
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Figure 1: Left: G\/Ap/Gpp vs. t for all masses which have been simulated. Right: (ap), obtained 
by fitting the data in the plot at left, as a function of the bare quark mass. The dashed line 
represents the result of a quadratic fit (see text). 



should not come as a surprise since, even if A^{x) is not the conserved current, it has the correct 
behavior under global non-singlet chiral transformations. 

In the quenched approximation the contribution of exact chiral zero modes of the Neuberger- 
Dirac operator is not suppressed by the fermionic determinant. In some correlation functions, this 
can give rise to large quenched artifacts for small masses. For example it is easy to show that 
Gpp{t) receives from the unsuppressed zero modes contributions proportional to and 1/m, 

which should vanish in the infinite volume limit, but can be quite sizeable for finite volume. A clever 



am 


Gs-p 


Gpp 




Zs-p 


aMp 


afp 


Zpp 


aMp 


afp 


0.100 


0.0040(4) 


0.379(6) 


0.089(2) 


0.0042(3) 


0.382(3) 


0.089(2) 


0.085 


0.0036(4) 


0.348(6) 


0.085(2) 


0.0039(3) 


0.352(4) 


0.085(2) 


0.070 


0.0033(4) 


0.315(7) 


0.081(2) 


0.0036(3) 


0.321(4) 


0.081(2) 


0.055 


0.0030(5) 


0.280(9) 


0.076(2) 


0.0034(3) 


0.287(5) 


0.077(2) 


0.040 


0.0026(5) 


0.239(11) 


0.071(2) 


0.0032(4) 


0.250(7) 


0.073(2) 



Table 1: Mesons masses and matrix elements for all the bare quark masses considered in the 
simulations, as obtained from Gpp{t) and Gs-p{t). 



way to avoid such artifacts has been proposed by the authors of Ref. [^], who noticed that the zero 
modes, because of their chirality properties, contribute equally to the Gpp and Gss correlation 
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Figure 2: Left: A representative example {am = 0.070) of the effective mass from Gs~p{t)- Right: 
as at left but from Gpp{t). 



functions, so that their contributions cancel in the difference 

Gs-p{t) = Gpp{t) - Gss{t) (33) 

which can also be used to extract the pseudoscalar meson mass and decay constant, since the 
contributions from the heavier scalar mesons fall off faster. The drawback is, of course, that the 
plateau in the effective mass and, correspondingly, the range that can be used for the cosh fit 
become shorter. Nevertheless we found the plateau in the Gs-p to be long enough to permit a 
meaningful fit. We fitted Gs-p{t) to a single particle propagator with a cosh dependence on t, 

Gs-p{t) = %^ exp(-^aMpr) cosh(aMp(5 - t)) (34) 
alVl p 2 2 

in the time interval ti — t2 = 12 — 16. The lower limit was fixed at the point where we found 
stabilization of the effective meson masses. The results of the fits are also given in Table || and 
an example of the effective meson mass from Gs-p{t) is shown in first plot of Fig. ^. We also 
performed a two cosh fit of Gs-p{t) finding consistent results. 

We have also fitted the correlation functions Gpp{t) to a single particle propagator with a cosh 
dependence on t, as in Eq. (^), in the time interval ti — t2 = W — 16. As before the lower limit is 
fixed as the point at which the values of the effective meson masses become stable. We report our 
results in Table |l] and in the second plot of Fig. ^ we give an example of the effective meson mass 
as extracted from Gpp{t). 

We illustrate in Fig. ^the values for (aMp)^, obtained from Gs~p{t) and Gpp{t), as a function 
of the bare quark mass am. In both cases a linear behavior 

{aMpf = Amp + BMp{am) (35) 
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Figure 3: (aMp)^ vs. am as obtained from Gs~p{t) (left) and Gpp{t) (right). The dashed hnes 
represent the results of a linear fit. 



fits very well the data with 



-0.0005(68) 



1.43(7) 



for the masses obtained from Gs-pit) and 



Amp = 0.006(4) Bmp = 1.39(3) 



(36) 



(37) 



for those obtained from Gpp{t). We see that the parameters of the two fits are compatible and that 
both intercepts vanish within statistical errors. In particular, one does not notice any sign of the 
singular contributions from zero modes in the masses obtained from the pseudoscalar correlation 
functions. These are expected to show up at some point, but one would probably need much higher 
statistical accuracy and lower values of m to bring them into evidence. 

Contrary to the pseudoscalar masses, the results for the matrix elements, parameterized as in 



Eq. (34) by the factors Zs-p and Zpp for Gs-p{t) and Gpp{t), respectively, show more significant 
differences. In Fig. ^we reproduce the results we obtained for Zs~p and Zpp respectively. We see 
that, while Zs-p exhibits a linear behavior as function of am, the graph for Zpp shows a clear 
indication of curvature. Linear fits of the form 



Zi = Ai + Bi{am) 



(38) 



give 



As^p = 0.0016(7) 
App = 0.0025(4) 



Bs-p = 0.024(6) 
Bpp = 0.016(3) 



(39) 
(40) 
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Figure 4: Left: Zg-p vs. (am); the dashed hne represents the result of a hnear fit. Right: Zpp 
vs. {am); the dashed hne represents a fit of the form of Eq. (p]). 



Within the quite large statistical errors, the intercept and the slope obtained from Gs^pit) are 
still compatible with those derived from Gpp{t). However the central values are quite different. If 
we interpret the curvature in the graph of Zpp as due to pole terms from the zero modes and try 
a simple fit of the form 

Vpp 

Zpp = App + Bppiam) + (41) 

am 

we obtain 

App = 0.0014(3) Bpp = 0.024(2) Vpp = 0.000035(16) (42) 

Now the x^/d.o.f of the fit turns out to be much smaller and the central values of intercepts and 
slope parameters for the fits of Zg-p and Zpp are much closer. The fact that the curvature shows 
up only in the results for Zpp points to the fact that what we are seeing is the effect of the 
unsuppressed zero modes, and not of chiral logarithms which would affect both sets of results (and 
would most likely become noticeable at much smaller values of am). On account of the above, we 
will use Gs-p{t) to derive our further results. It must also be said that some of the observables will 
be calculated directly at m ~ ms/2 (see below), and for these observables the difference between 
Zpp and Zs-p is irrelevant within our statistical error. 



The pseudoscalar decay constant is defined through Eq. ( |22D and we calculated its value directly 
from the parameters extracted from a cosh fit to Gs-pit) 



The values we obtained for afp in correspondence to each simulated mass are reproduced in Table |^ 
(where we also reproduce the values derived from Gpp{t)) and are shown in Fig. |5|. The dashed and 
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Figure 5: The procedure used to derive afx and aMx- The dashed and dashed-dotted hnes 
represent fits to our numerical results. The solid line represents the curve (a/p) = Csi{aMp), 
where the coefficient Csi = fj^^/M^^. 



the dot-dashed lines show the results of a linear and a quadratic fit in (aMp)^, respectively. The 
data slightly favors a quadratic fit, but our statistical accuracy is insufficient to rule out a linear 
fit as inconsistent with the data. It is interesting to observe, though, that we obtain fx/f-K — 1.14 
and ~ 1.23 from the linear and quadratic fits, respectively, with the results from the quadratic fit 
in much better agreement with the experimental value. 

Starting from our values for Mp and fp, we fix the lattice spacing and the physical meson 
masses by using the method of "lattice physical planes" |^]. This avoids recourse to a chiral 
extrapolation for observables where it is not really needed. In the plane [a/p, (aMp)^], we plot our 
lattice data as well as the curve (a/p) = Csi{aMp) (the solid line in Fig. ||), where the coefficient of 
proportionality Csi {si for strange-light) is fixed by the experimental value for the ratio /M^^ . 
^ The point where the two lines meet determines afx as well as aMx- Since our statistical accuracy 
does not allow us to discriminate between a linear and quadratic fit to the data for afp vs. {aMp)"^ 
, we just use the linear fit for the tiny extrapolation to the "Kaon" point (see Fig. |5|). Our results 
are 

aMx = 0.216(8) afx = 0.0698(26) (44) 
(A quadratic fit produces a result which differs from the above by far less than the statistical 

•^Throughout the paper we use the following experimental numbers: M'^^ = 495 MeV, f1^^ — 160 MeV. We 
neglect the experimental errors which are well below our statistical and systematic errors. 
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errors). To set the lattice spacing we compare the value of afx with its experimental value and we 
get 

aj^ = 2.29(9) (45) 

This is the value of the lattice spacing we will use throughout the paper. Again we stress that up to 
this stage we did not have to perform any chiral extrapolation, but derived all physical information 
staying close to the region of quark masses used in the actual calculations. 
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Figure 6: The data used to derive the chiral condensate from the GMOR relation. 



From our data, we can also compute the chiral condensate by using the GMOR relation of 
Eq. (p]). In Fig. ^ we show the quantity 



-'^rn = -{am)^^^ (46) 



as a function of the bare quark mass. The data exhibits a very good linear behavior and a linear 
fit leads to 

a^;)^ = -0.00117(27) (47) 

We performed also a quadratic fit to the data, obtaining a result consistent with zero for the 
coefficient of the quadratic term, but with an error twice as large in the intercept, because of the 
additional degree of freedom. Therefore we conclude that we see no indication for a quadratic term 
in the fit within our statistical error. 



The chiral condensate can also be computed directly from Eq. ([Tq). However, away from the 
chiral limit, the r.h.s. of that equation contains power divergent terms (in a~^, see Eq. (pO|)) and, in 
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the quenched approximation, infrared divergent contributions (in m) due to the unsuppressed zero 
modes. Therefore quite a severe extrapolation from our simulated data is needed to determine its 
chiral value. One can take advantage, though, of a computational strategy similar to that used in 
Eq. (|33D. To remove the contribution due to unsuppressed zero modes, we subtract the scalar-scalar 
correlation function from the r.h.s of Eq. (jl^). From a quadratic fit in the quark mass we obtain 

a^xia) = -0.00117(42) (48) 

which is in remarkable agreement with Eq. (^7[). Of course, the larger error is an indication of 
the difficulty of the extrapolation. Larger statistics and a careful treatment of the zero mode 
contributions would be needed to obtain more reliable and precise results by the direct method. 



4 Physical Results 

In this section we will use our lattice results to infer the renormalized values of the sum of the 
strange and average up-down quark masses and of the chiral condensate. To obtain our final 
values, we still need the value of the the scalar renormalization constant Zs{fJ^a), which we have 



computed non-perturbatively in the RI/MOM scheme following the approach proposed in |17|. The 
details of these calculations will be presented in a separate paper |2^]. The result we obtained in 
the RI scheme is 

Zfi(2GeV) = 1.24(5) (49) 

where the error is mainly systematics due to the chiral extrapolation and to the uncertainty in 
the value of the renormalization scale^ Judging from the comparison of the data at different 
renormalization scales with the logarithmic evolution predicted by the renormalization group equa- 
tions at N^LO, the discretization errors appear to be well below the error in Eq. (|49|). The result 
in Eq. ( ^9|) differs by about 10% from the bare perturbation theory determination at one loop 
\Zg-^{2GeV)]^'^ = 1.11 [^, 24 1 . However if we use as expansion parameter we obtain a result 



consistent within errors with Eq. (^91). By using N^LO continuum perturbation theory |^3[ with 
Nf = and Aqcd = 0.238(19) ||ll| to convert the result in Eq. (^) into the MS scheme, we obtain 



Z|^S(2GeV) = 1.41(6) (50) 

Had we used the experimental N^LO results from as{Mz) = 0.118, we would have obtained a value 
of the scalar renormalization constant ~ 10% higher than the one given above. The difference can be 
taken as an indication of the systematic error introduced by the quenched approximation. Therefore 
we will add this uncertainty in quadrature to our final estimate of the systematic error. N^LO 
perturbative computations are available for the RI — MS matching of the scalar renormalization 



constant |34], but the difference with the N LO results is much below the error induced by the 



quenching ambiguity discussed above. If we had used the procedure proposed recently in Ref. |27] 



on our data, we would have obtained a value for Zg^^^(2GeV) which agrees within the errors with 
the one in Eq. (|50|), but with a central value ~ 10% higher. 



We stress again that this number is obtained at /9 = 6.0 for p — 1.4 and r — 1.0. 
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We finally combine this last result with the numbers presented in the previous section to 
obtain our main physical results. From the value of aMx in Eq. ( p^ ) and from the lattice spacing 
in Eq. (|4^) , we obtain for the combination of bare quark masses 

ms{a) + m(o) = 149(9) MeV (51) 

where the error is only statistical. Since our volume is fairly large, we expect our main sources of 
systematic errors to come from discretization effects and from the quenched approximation. For 
a rough estimate of the systematic error due to quenching approximation we can use the results 



in Ref. |36|. From these one sees that, within the quenched approximation, the use of different 
observables to calculate the lattice spacing can produce differences in the results of ~ 10%. Had 
we used tq to fix our lattice spacing we would have obtained a number ~ 7% higher than the 



one in Eq. (pl|). Combining this fact with the results in |12|, we would infer that discretization 
uncertainties are below 10%. In order to be conservative, we will take 15% as the estimate of 
our overall systematic error in the renormalized quark masses due to quenching and discretization 
effects. A more precise estimate of the systematic errors will need much more extensive simulations, 
which at present would be beyond our capability and the exploratory scope of this work. Combining 



the results in Eqs. (50) and (51) we obtain 

(m, + m)^^(2 GeV) = 120±7±21MeV (52) 

which represents one of the major results of this paper. By using Eqs. @ and §^ the above 

translates to 

mfS(2 GeV) = 102±6± 18 MeV (53) 



This result agrees very well with the current lattice world average [16|. 



Insofar as the value of the condensate is concerned, if we used the standard two-step approach. 



i.e. first measure the dimensionless condensate, see Eqs. (47) or (p8|), and then multiply it by the 
cubic power of the lattice spacing, the result would be affected by a very large systematic error due 
to the uncertainty in the determination of the lattice spacing in quenched simulations. Instead, 
we will use an alternative method Q. We write the GMOR relation (|2^ ) for the renormalized 
condensate as follows 

X(a) = -l/2SMpa-i (54) 



where Bmp is defined in Eq. (|36f ) and = 0.1282 GeV is the "experimental" value, in physical 
units, of the pseudoscalar decay constant extrapolated to the chiral limit. While computing the 
condensate from the above formula relies on an additional element of experimental information, it 
has several advantages. The most important is that, by expressing the condensate in terms of f^, 
we are left with only one power of the UV cutoff a~^. Another advantage is that, if we assume that 
the relation between Mp and m stays substantially linear for a range of quark masses extending to 
~ TTT-s (and our numerical results validate this, see Fig. ^), then there is no need for an extrapolation 
to the chiral limit to evaluate Bmp- With this method we obtain 

{■ipil))^^ (2 GeV) = - 0.0167 ± 0.0010 ± 0.0029 GeV^ (55) 
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where the estimate of the systematic error has been made using the same criteria we used for 
estimate of the error in the quark mass. This is our best value for the chiral condensate. It is 
interesting to note that, if we had used the standard technique, starting from Eq. ( p7|) and with 
the lattice spacing of Eq. (pSj), we would have obtained a result with a central value very close to 
the value in Eq. (|55|), but with a much higher systematic error. 

Finally, from Eq. (p5|) and N^LO matching, we get 

{i'^pf^{2 GeV) = - 0.0190 ± 0.0011 ± 0.0033 GeV^ = - [267 ± 5 ± 15 MeV]^ (56) 
This result is in very good agreement with the result obtained by the authors of Refs. []35|, [27[|, 



while it is smaller than the result in Ref. [38|, even if still compatible within errors. Our result is 
also compatible within errors with the number obtained few years ago in Ref. Q with Wilson-type 
fermions. We expect, though, the systematics due to the discretization effects to be smaller (0(a^)) 
in the result reported in Eq. (56) than the error {0{a)) which affects the determination in Ref. 0]. 



5 Conclusions 

In the overlap regularization chiral symmetry is preserved at finite lattice spacing and finite volume, 
therefore there is no mixing among operators of different chirality. As a consequence no additive 
quark mass renormalization is required and no fine tuned parameters are needed to compute 0{a) 
improved masses and matrix elements. Our results have indeed produced a remarkable verification 
of "good chiral behavior" both in the axial Ward identity and in the pseudoscalar masses. 

In this paper we presented the results of the first computation of (ms+m) with overlap fermions 
in the quenched approximation. To avoid uncertainties due to lattice perturbation theory, we 
computed the multiplicative renormalization constant ZsifJ^a) non-perturbatively in the RI/MOM 
scheme. Our main results have been summarized in the introduction. While the systematics 
errors due to quenching are common to previous calculations, the other systematic errors (mostly 
discretization effects) are different than with other lattice regularizations and likely to be smaller, 
because of chiral symmetry. We also computed the chiral condensate (i/^ip) from the GMOR relation 
and directly. Even if, given our statistical and systematic errors, the former method is more reliable, 
it is rewarding to notice that the two determinations are in good agreement. 

The calculation of light quark masses uses many of the ingredients needed for a lattice calcula- 
tion of weak matrix elements, although the latter is computationally more demanding. From this 
point of view, the very good agreement between our results for the quark masses and the current 
lattice world average bodes well for the use of the overlap formalism also for matrix element calcu- 
lations. Our investigation has been mostly of exploratory nature. One would need to extend it to 
larger volumes and better statistics. One should also find a more direct way to isolate and account 
for the effects of the zero modes. Nevertheless, we believe that it gave a strong indication that 
the overlap formalism can be used effectively, with known algorithms and the present generation 
of computers, for large scale QCD calculations, at least in the quenched approximation. Thus we 
would conclude that it represents a very promising non-perturbative regularization for solving long 
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standing problems, such as the proof of the A/ =1/2 rule and the calculation of e'/e, which would 
be hard to address with conventional regularizations. 
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